Stable First-order Particle-frame Relativistic Hydrodynamics for Dissipative Systems 
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We propose a stable first-order relativistic dissipative hydrodynamic equation in the particle frame 
(Eckart frame) for the first time. The equation to be proposed was in fact previously derived by 
the authors and a collaborator from the relativistic Boltzmann equation. We demonstrate that 
the equilibrium state is stable with respect to the time evolution described by our hydrodynamic 
equation in the particle frame. Our equation may be a proper starting point for constructing 
second-order causal relativistic hydrodynamics, to replace Eckart's particle-flow theory. 
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Relativistic hydrodynamics (RHD) is a useful tool for 
analyzing slow and long wavelength behavior of rela- 
tivistic many-particle systems in terms of static and dy- 
namic thermodynamic properties. In fact, RHD is widely 
used in astrophysics [1] and the phenomenology of rel- 
ativistic heavy ion collisions Since works demon- 
strating the success of perfect hydrodynamics in de- 
scribing the phenomenology of the Relativistic Heavy 
Ion Collider (RHIC) at BNL 0, i, H, we are witness- 
ing agrowing interest in RHD for dissipative systems 
[IBUIll- Indeed, there have been many works attempt- 
ing to show how small can be the transport coefficients 
of strongly-interacting systems composed of hadrons or 
quarks and gluons, with many of these employing the 
so-called AdS/CFT correspondence hypothesis f9'|. It 
should be noticed, however, that the theory of RHD for 
dissipative systems is not clearly established, although 
there have been many fundamental studies since Eckart's 
pioneering work [Toj . 

We identify the following three fundamental problems 
regarding relativistic hydrodynamic equations (RHDEs) 
for dissipative fluids I 111 : (a ) ambiguities in the forms of 
the equations [1, 0, Efl [Hj [Hj [I3l; (b) the unphysical 
instability of the equilibrium state in the theory of the 
so-called first-order equations, in particular in the Eckart 
frame [l^ , defined below; (c) the lack of causality in the 
first-order equations [H, [H, [13, [l3|. The present pa- 
per is concerned with the first two problems. Although 
the unphysical instability of the equilibrium state may 
be attributable to the lack of causality, and the Israel- 
Stewart equations with second-order time-derivative are 
presently being examined in connection to this problem 
p, H, 0, [ill, we emphasize that the first two problems 
and the third one have different origins, and the first 
two must be resolved before the third is addressed. Note 
that the causality problem also exists in non-relativistic 
cases and is in essence a problem of how to incorporate 
the space-time scales shorter than those corresponding to 
the mean-free path, beyond those in the usual hydrody- 
namic regime. We also remark that the proper form of 
Israel-Stewart-type equations has not yet been definitely 
determined @, |3|- 

Let us represent the flow velocity by u^, with — 



gi"'u^u^ = 1 {gi"' = diag(-hl, -1,-1,-1)). In the rel- 
ativistic theory, the rest frame of the fluid and the flow 
velocity cannot be uniquely deflned when there exist 
viscosity and heat conduction. In the phenomenological 
theories [l^, [i2|j the ambiguity of the flow velocity 
is resolved by placing constraints on the dissipative part 
of the energy-momentum tensor, 5T^'^ , and the particle 
current, i5iV^. Landau and Lifshitz deflned such that 
there is no dissipative energy density, energy flow nor 
particle density; i.e., we have the constraints (JT^" = 
(referred to as ET) and -u^ SNf" = (EN). This frame is 
called the energy frame. Contrastingly, Eckart chose the 
particle frame, in which there is no dissipative contribu- 
tion to the particle current; i.e., we have 5N^^ — (PN), 
together with Ui, ST^''' = (PT): These conditions im- 
ply that there is no dissipative contribution to the energy 
density in this frame. However, it should be noted that 
the seemingly plausible constraint PT on ST^'^ is prob- 
lematic, as shown in 11] and explained below. 

Recently, Tsumura, Kunihiro (the present authors) 
and Ohnishi (abbreviated as TKO) [11[ derived generic 
covariant hydrodynamic equations for a viscous fluid 
through a reduction of the dynamics described by the 
relativistic Boltzmann equation in a systematic manner, 
with no heuristic arguments, on the basis of the so-called 
renormalization group (RG) method HO, [H. This 
was done by introducing the macroscopic frame vector 
that defines the macroscopic Lorenz frame, in which 
the slow dynamics are described. The generic equation 
derived by TKO can produce a relativistic dissipative hy- 
drodynamic equation in any frame with the appropriate 
choice of a^; the resulting equation in the energy frame 
coincides with that of Landau and Lifshitz [l2|, while 
that in the particle frame is similar to, but slightly dif- 
ferent from, the Eckart equation. Interestingly, the TKO 
equation in the particle frame does not satisfy the con- 
straints PT on ST^"" but, instead, satisfies STt'^ = 0, 
which we call PT', together with PN. It should be noted 
that the new constraints, PT', are identical to a match- 
ing condition postulated by Marie and Stewart (MS) in 
the derivation of the RHD from the Boltzmann equation 
with use of Grad's moment theory ^| . We call the con- 
straints PT', together with PN, the Grad-Marle-Stewart 
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(GMS) constraints. In [Tl|, TKO proved that the simul- 
taneous constraints PT and PN cannot be compatible 
with the underlying Boltzmann equation if the hydrody- 
namic equation describes the slow, long wavelength limit 
of the solutions of the Boltzmann equation. This is in- 
teresting in connection to problem (b), i.e., the fact that 
the solutions of the Eckart equation around the thermal 
equilibrium are unstable , while the Landau theory is 
stable. 

An immediate question is whether the solutions of the 
new equations in the particle frame are stable around the 
thermal equilibrium. In fact, the hydrodynamic equa- 
tions of MS and TKO in the particle frame are of differ- 
ent forms, although both satisfy the constraints PT' and 
PN. In the present paper, we examine the stability prob- 
lem for the new equations in the particle frame. Because 
second-order equations, such as the Israel-Stewart equa- 
tions, are usually constructed in the particle frame, as an 
extension of the Eckart equation, finding a stable first- 
order equation in the particle frame is of fundamental sig- 
nificance. As the RG method has been employed to con- 
struct the slow dynamics of various systems through the 
explicit construction of the slow, stable manifold of the 
dynamics, we conjecture that the hydrodynamic equa- 
tion obtained as the slow, long wavelength limit of the 
Boltzmann equation on the basis of the RG method will 
provide a description in which the thermal equilibrium 
state is stable. We demonstrate that this is indeed the 
case by performing a linear stability analysis using the 
EOS and the transport coefficients for a rarefied gas. By 
contrast, we find that the MS equation, like the Eckart 
equation, is unstable. Hence, for the first time, a stable 
RHDE is obtained in the particle frame. We believe that 
this will provide a sound starting point for the construc- 
tion of the proper second-order equations. 

The energy-momentum tensor for our equation in the 
particle frame reads 

+ C{?,u^'u'' - /^^"') [-(37- 4)-2 V • u] 
-H77(V^u''-f V''u''~2/3A^''V-u), (1) 



while the particle current is given by N^^ 



i^, with 



9 



u'^u" and V = A^-^a^. Here T, ^, e, 
p, n and 7 are the temperature, the chemical potential, 
the internal energy, the pressure, the particle density and 
the ratio of the specific heats, respectively, and C, A and 
•q denote the bulk viscosity, the heat conductivity and 
the shear viscosity, respectively. The MS equations are 
obtained from the above equations through the replace- 
ments -C(3 7 - 4)"2 V • u — > +C(3 7 - 4)"i V • u and 
AV^T — > A(V^T - TDut'), where D = u" d^. One 
can easily check that both equations satisfy the GMS 
constraints. Nevertheless, we find the following differ- 
ences between them: (A) the thermal forces in the MS 
equations contain the time-like derivative of the flow ve- 
locity Du^ , while those in our equations involve only the 
space- like derivative V, and (B) the sign of the thermo- 
dynamic force owing to the bulk viscosity in our equation 



is the same as that in the Landau equation and opposite 
that in the MS equation. We can trace the two charac- 
teristic features of our theory back to the simple ansatz 
that only the spatial inhomogeneity, over distances of the 
order of the mean free path, is the origin of the dissipa- 
tion. It should be noted that the same ansatz for the 
non-relativistic case leads naturally to the Navier-Stokes 
equation, as shown in [2l|, and hence our framework can 
be interpreted as the most natural covariantization of the 
non-relativistic case. 

The thermal equilibrium state is given by u'^{x) ~ 
(1, 0, 0, 0) = Mq, T{x) = To and — /io, with Tg 

and /io being constant. This is a trivial solution to 
the equations. Let us investigate the linear stability of 
the equilibrium solution. Writing T(x) = Tq + ST{x), 
li(x) = Ho + S^j,{x) and u^^{x) = Uq + Su^{x), we ex- 
amine the time evolution of the deviations in the lin- 
ear approximation using the evolution equation given by 
a^T'"' = and d/.N'^ = 0. Here we note that the in- 
dependent variables are the five quantities 6T[x), 5ijl{x) 
and 5u'^{x) {i = 1, 2, 3), because 5vP{x) = 0, due to the 
constraint u^{x) = 1. 

In terms of the Fourier components $Q(fc) = 
*(5wi(fc), Su'^ik), 6u^{k), Sfik), Stl{k)), defined through 
^a{x) = /^l^ $Q(fc) e"**^'^, the linearized hydro- 
dynamic equation reduces to the algebraic equation 
ELi Mc^P = 0, with 
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where we have set A:^ = (fc", 0, 0, fc'^) without loss of 
generality. The first and second components of $q. de- 
scribe the transverse mode, while the third component 
the longitudinal one. Here Ci^i^io are given by Ci = 
{e + p) i-i fc") + 7/ |fe|2, £2 = -?//3 - Cp, = dp/dT - 
A(-zfcO), £4 = dp/dfi, £5 = -ie+p)+3Cpi-ik°) 
£6 = de/dT i^i fcO) + A \k\^, £7 = de/d^i {-i k°), £§ = 
— n, £9 = dn/dT {—i k°) and £10 = dn/dfi{—i k°), with 
(^P = C (3 7 — 4)"^ being the effective bulk viscosity in the 
particle frame. In the above, the quantities, e, p, n, 7, ^, 
A, ry, de/dT, de/dp, dp/dT, dp/d\i, dn/dT and dn/dpL 
take their equilibrium values, with T = Tq and jjL = pLQ. 

The existence condition of a solution reads det M = 0, 
which reduces to 

£? [(£1 - |fcp £2) (£6 £10 - £7 £9) - Ifcp £5 (£3 £10 

-£4£9)-|fcp£8(£4£6-£3£7)] =0. (3) 

This equation gives the dispersion relation fc° = fe'^dfej) 
for the hydrodynamic modes, and the stability condition 
for the equilibrium state reads Imfc° < 0, V |fc|. 

We see the dispersion relation for the transverse mode 
is given by £1 = 0, whose solution is k^ = —ir] |fcp/(e -I- 
p). Thus, we find that the transverse mode is stable. 
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Here we again stress that the equation we study does 
not contain a term proportional to Du'^ in the thermal 
force for the heat flow. What would happen if such a 
term were present in the thermal forces, as in the case 
of the MS and the Eckart theories? In this case, the 
corresponding equation becomes £i = (e + p) {—i — 
TX{—ik°)^ +7y|fcp = 0, which possesses a root with 
a positive imaginary part, and hence an unstable trans- 
verse mode appears. We emphasize that this instability 
is inevitable if the heat flow term contains Du'^ 0]. 

Next, we examine the dispersion relations of the lon- 
gitudinal modes. We first consider the simple but inter- 
esting case in which the heat conductivity vanishes (i.e., 
A = 0), but the bulk and the shear viscosities may be 
positive (i.e., C 7^ ^nd r] ^ 0). This simple case is 
often studied in the literature. We subsequently carry 
out a full analysis in which all the transport coefficients, 
including A, may be positive. 

In the simple case with A = 0, the equation has the root 

— and those satisfying oq {—i k^)^+bo (— i A;") + co = 

0, with ao = (e+p) {e , n}, bo = [(4r//3 + Cp) {e , n}- 
^Cp{P,n'}] and Co = \k\'^ [{e + p) {p , n} + n {e , p}], 
where we have written the Jacobian as {F , G} = 
d{F , G)/d{T , /i). Now, a simple analysis of the alge- 
braic equation aouj'^ + bauj + cq — with uj — ^ik'^ 
shows that the necessary and sufficient condition for 3 k'^ 
with Imfc" < is that ba/ao > and cq/qo > 0. Owing 
to the properties of the Jacobian and the thermodynamic 
relations, the last inequality generally holds, because the 

1. h.s can be rewritten as co/ag — \k\'^ {dp/de)s — |fe|^Cg, 
with Cs being the sound velocity. Then, the stability con- 
dition reduces to bo/ao > 0, which is equivalent to 

Ar,/3 + Cp [1 - 3 {dp/deU > 0. (4) 

This is a new condition that involves not only the EOS 
but also the bulk and shear viscosities. It can be shown 
analytically [2^ that this inequality is satisfied at least 
for a rarefied gas in the massless limit. To see this, first 
notice that e = 3 p for a relativistic gas composed of mass- 
less particles. Then the inequality reduces to the trivial 
one rj > 0, because the second term with a bracket on 
the l.h.s vanishes, provided that the effective bulk vis- 
cosity, Cp = C (37 — 4)"^, is finite in the massless limit. 
In fact, it can be shown that this is the case using the 
microscopic formula for ( [2^ , although 3 7 — 4 ^ in 
the massless limit. We also remark that numerical cal- 
culations using the viscosities C and 77 obtained from the 
Boltzmann equation reveal that the inequality ([4]) is al- 
ways satisfied, even for a rarefied gas of massive particles. 
Instead of presenting the numerical results for this limit- 
ing case, we present the results for the general case, i.e., 
that in which A 7^ 0, C 7^ and ry 7^ 0, below. 

We now demonstrate that the thermal equilibrium 
state is stable with respect to the dynamics described 
by our equation, even when the heat conductivity. A, 



is finite. The dispersion equation for the longitudinal 
modes is obtained from the roots of the cubic equation 
a Lj'^ + b uj'^ + c Lj + d = 0, with uj = —i /c°, where the coeffi- 
cients are given by a = ao-l-|fcp 3 Cp A {dn/dfi)T, b = 60 + 
nX\k\^ {de/dn)T, c = cq + \k\^ (47//3 + Cp) A {dn/dfi)T 
and d = n A (9p/9/i)T- The condition Imfc° < im- 
plies that the above equation for w has roots only in the 
left half plane or on the imaginary axis in the complex lo 
plane. An elementary analysis shows that this condition 
is given by 

a > 0, 5 > 0, d > and 6c- ad > 0. (5) 

Here, the equality holds in the case that the imaginary 
part of fc" vanishes. Note that the above equalities imply 
that c > 0. 

Now we demonstrate that these inequalities are 
satisfied for rarefied systems. For a relativistic 
free gas, we have n = {2^)'^ Amn^ ct [z~^ K2{z)], 
€ = mn[K^[z) / K2{z) — p = nT and 

7 = I + [z"^ + - {h - with z ^ m/T and 

h = {e+p)/nT being the reduced enthalpy. Here, K2{z) 
and K3{z) denote the second and third modified Bessel 
functions, respectively. It is seen that the positivity 
condition 0? > holds from the formula p = nT 
with n oc eTj which implies that {dp/dfj,)T > 0. It 
remains to demonstrate the rest of the inequalities, 
i.e., a > 0, 5 > and be — ad > 0, for which we 
need explicit forms of the transport coefficients as well. 
The transport coefficients Ci A and rj for a rarefied 
gas can be obtained from the collision term in the 
Boltzmann equation. The Galerkin approximation 
using the Ritz polynomial expansion fl7] with a con- 
stant cross section a in the collision integral gives 
C = m^^^''^^[^^'Kl{z)[{5 - 3j)h - 37]2]/[2if2(2z) + 
zK,{2z)], A = J-ieM/T[^2^2(^)[^/(^ _ 

l)]2]/[(z2 + 2)K2{2z) + 5zK3{2z)] and r, 
^^e''/^[z2i^|(z)/i2]/[(5z2 ^ 2)K2{2z) + {3z^ + 
4:9z)K3{2z)]. Note that all the transport coefficients 
are proportional to the inverse of the cross section, a. 
This implies that a strongly (weakly) interacting system 
has small (large) transport coefficients. The numerical 
results for a, b and be — ad are displayed in Figll] 
where the z = m/T dependence is shown using aT^ = 1 
for a wide range of values of the three momentum: 
\k\/T = 0.1 - 10. We have confirmed that the positivity 
of these quantities holds for a wide range of values of 
the cross section: aT^ = 0.01 - 10. We point out that 
a rarefied gas is a system in which dissipative effects 
are most significant. Thus, we have demonstrated that 
the thermal equilibrium solution is stable within the 
description provided by our hydrodynamic equation 
in the particle frame. Obviously, a solution with flow 
is unstable in a viscous fluid, as it must relax to the 
equilibrium state. 



FIG. 1: The 2 = m/T dependence of a/(r^e^''/'^), h/{T^° e'"''^) and {he - ad) / T^) for \k\/T = G.l, 1, 10. Some factors 
have been multiphed so that the variables become independent of /i. It is seen that all these quantities are positive in all cases. 



In this Letter, we first pointed out that the constraint 
proposed by Eckart, u^Ui/dT^^" = (PT), is incompati- 
ble with the underlying relativistic Boltzmann equation. 
This important point has been largely unnoticed. We 
then showed that the reduction of the Boltzmann equa- 
tion employing the RG method leads to an RHDE in 
the particle frame, that satisfies the constraint 5T^^^ — 0, 
while Ui, T'"' = e — 3 V • u, which includes a contri- 
bution from the flow as well as the internal energy. This 
equation might imply that the energy density of an ex- 
panding system extracted from the hydrodynamic anal- 
ysis can be erroneous. We have demonstrated that the 
solution around the equilibrium state in the new equa- 
tion is stable. This was done by carrying out a linear 
stability analysis using the EOS and the transport co- 



efficients for a rarefied gas. We conclude that Eq. (1) 
represents the first viable possibility as a stable, first- 
order, particle-frame RHDE for a viscous fluid. This is 
significant because the Israel-Stewart causal equation is 
usually constructed in the particle frame with PT. A de- 
tailed presentation of this work and applications of the 
new equation studied here will be reported elsewhere. 
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